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We investigate the response behavior of an Ising system, driven by an oscillating field, on a small- 
world network, with particular attention to the effects of the system size. The responses of the 
magnetization to the driving field are probed by means of Monte Carlo dynamic simulations with 
the rewiring probability varied. It is found that at low and high temperatures the occupancy ratio, 
measuring how many spins follow the driving field, behaves monotonically with the system size. At 
intermediate temperatures, on the other hand, the occupancy ratio first grows and then reduces as 
the size is increased, displaying a resonance-like peak at a finite value of the system size. In all 
cases, further increase of the size eventually leads to saturation to finite values; the size at which 
saturation emerges is observed to depend on the temperature, similarly to the correlation length of 
the system. 
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It has been known that under appropriate circum- 
stances, the presence of noise in a driven system can 
enhance rather than suppress the output of the system. 
Such attractive phenomena, called stochastic resonance 
(SR), have been widely investigated in various systems 
of many practical applications P, @] . Recently, those SR 
phenomena have also been examined in the context of 
the system size resonance: Stochastic flips of the mean 
field are observed to depend on the system size, leading 
the linear response of the system to reach a maximum 
at a certain system size Q. Those works considered two 
systems of fully coupled noisy oscillators and one two- 
dimensional Ising system with nearest neighbor interac- 
tions. Namely, the underlying connection topology of dy- 
namic variables was assumed to be regular, with either 
local or global connections. Meanwhile, recent studies of 
neuronal networks, computer networks, biochemical net- 
works, and even social networks, have revealed that var- 
ious real systems in nature possess quite complex struc- 
tures, which can be described neither by regular networks 
nor by completely random networks iJj . Accordingly, it is 
desirable to study effects of the system size on collective 
responses in the systems with the connection topology 
of complex networks, which can be made more realis- 
tic. In particular, the interplay between the system size 
and noise may be relevant in various biological systems 
such as neural networks and other cell networks, which 
consist of finite numbers of elements. For example, in 
the study of stochastic resonance in biological systems, 
optimal sizes of calcium ion channel clusters have been 
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examined. Observed is that the clustering of the release 
channels in small clusters increases the sensitivity of the 
calcium response '.5\ . This suggests a possible realization 
of the system, providing motivation for the investigation 
of the size resonance in the complex network structure. 

In this paper we consider an Ising model on a complex 
network, specifically, on the Watts and Strogatz (WS) 
type small- world network |^. It is well known that the 
WS network is characterized by a small characteristic 
path length £ ~ In TV, where N is the number of nodes 
constituting the network, and a large clustering coeffi- 
cient. Both are commonly observed properties of real 
networks in nature. The WS network in this paper is 
constructed following Ref . |^ : A regular one-dimensional 
(ID) network of N nodes is first constructed with local 
connections of range k, under periodic boundary condi- 
tions. At this stage each node on the network has 2k 
nearest neighbors. Next, each local link is visited once 
and, with the rewiring probability P, removed and recon- 
nected to a randomly chosen node. After a whole sweep 
of the entire network, the average number of shortcuts in 
the network of size N is given by NPk. Accordingly, the 
rewiring probability P may be regarded as the fraction 
of the average number of shortcuts over the total number 
of connections Nk. In this paper, the local interaction 
range k is set equal to two for convenience; longer ranges 
{k > 2) are not expected to lead to any qualitative dif- 
ference. After the WS network is built as above, an Ising 
spin is put on every node, and an edge (or a link) con- 
necting two nodes is regarded as the coupling between the 
two spins at the two nodes. Finally, we apply an oscilla- 
tory field, driving the Ising spins, and the corresponding 
responses of the average spin, i.e., the magnetization are 
probed via Monte Carlo dynamic simulations, with at- 
tention to the effects of the system size. 
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The Hamiltonian for the field-driven Ising model on 
the WS network, which is constructed as described above, 
reads 

H = -^^Jijaiaj - h{t)^a^, (1) 

where the ferromagnetic spin-spin interaction strength 
Jij is given by 

r - 7 ^ / J for j e Ai, , , 

"^'^ ~ - \ otherwise. ^ ' 

The neighborhood of node i stands for the set of nodes 
connected to i (via either local edges or shortcuts), and 
(Tj (= ±1) represents the Ising spin at node i. The si- 
nusoidally oscillating magnetic field hit) = ho cos ilt is 
applied with the driving amplitude ho and frequency Q, 
while the system is assumed to be in contact with an 
isothermal heat bath at temperature T. We probe the 
dynamics of the system described by Eq. ^ by means of 
Monte Carlo (MC) dynamic simulations, employing the 
heat bath single spin-flip algorithm and measuring the 
time t in units of the MC time step. For thermalization, 
we start from sufficiently high temperatures and lower 
the temperature T slowly with the decrement AT = 0.01 
(in units of J/ks with the Boltzmann constant ks)- The 
driving amplitude and frequency are taken to be ho = 0.1 
and = 0.001. We have also considered different fre- 
quencies, for example, = 0.01 and 0.1, and found that 
the resonance-like peak indicating the system size reso- 
nance behavior tends to diminish at higher frequencies 
(see below). While simulations are performed at a given 
temperature, the data from the first 10^ MC steps are 
discarded, which turns out to be sufficient for stationar- 
ity, and measurements are made for next 10^ MC steps. 
Networks of various sizes are constructed as described 
above, and averages are taken over 100 different network 
realizations for each size. 

To investigate the collective response of the system, 
we measure the occupancy ratio R which is defined to be 
the average fraction of the spins in the direction of the 
external field 8, 9J: 

^ _ / number of spins in the direction of h{t) \ 
\ total number of spins / ' 

where (• • • ) denotes the time average. In other words, 
R measures how many spins follow the oscillating mag- 
netic field. It is easy to understand that R approaches 
the value 1/2 in both low- and high-temperature limits 
(see, e.g., Ref. H [l3) and becomes increased near the 
stochastic resonance temperature, reflecting that more 
spins follow the external driving. Such SR phenomena 
have been observed in the system of given size, and it 
has been demonstrated that the matching condition of 
two time scales, the relaxation time of the system and 
the inverse frequency of the driving field, yields the op- 
timal noise strength Tsr at which the system displays 



maximum responses [9l llO|. Here we consider the sys- 
tem studied in Ref. [l^l from a different point of view 
and examine the behavior of the occupancy ratio with 
the system size at various temperatures, probing the size 
resonance. 

In the absence of long-range shortcuts (P = 0), the 
network structure reduces to that of the ID regular net- 
work with only local couplings of range k. Accordingly, 
when P = 0, the field-driven Ising model described by 
Eq. as well as the undriven model [h{t) = 0] does not 
exhibit long-range order at finite temperatures, yielding 
Tc = 0. Note, however, that even such a ID system 
displays SR behavior at finite temperatures j^J. In the 
presence of long-range shortcuts (P ^ 0), on the other 
hand, it has been found that both driven and undriven 
Ising model display ferromagnetic order at finite tem- 
peratures 0,01; furthermore, double SR peaks, which 
originate from matching of two time scales have been ob- 
served. 

We first consider the case without long-range shortcuts 
(P = 0), i.e., the purely ID system, and show in Fig.nthe 
behavior of the occupancy ratio R versus the system size 
N at various temperatures. It is observed that the occu- 
pancy ratio R first increases and eventually saturates to 
a finite value as the system size N is increased. The satu- 
ration size Ns , beyond which the occupancy ratio R does 
not increase any more, reduces as the temperature T is 
raised. Figure |2] displays such temperature dependence 
of the saturation size Ng, which has been taken as the 
size giving the occupancy ratio with the difference from 
the stationary value less than 5 x 10^**. One can observe 
the exponential behavior: Ns oc e"^/"^ with c = 5.6 ± 0.7, 
which is reminiscent of the behavior of the correlation 
length. In the ID Ising model the correlation length £ di - 
verges in the low-temperature limit as £ ~ gk(k+i)/T 
Here the local interaction range k has been chosen to be 
two (fc = 2), leading to the behavior £ ^ e^/-^, essentially 
the same as that of the saturation size. It is thus con- 
cluded that the saturation behavior emerges when the 
system size reaches the correlation length of the system. 

Meanwhile, in the presence of long-range shortcuts 
(P 7^ 0), substantially different behavior has been ob- 
tained for the occupancy ratio. In Fig. |31 the occu- 
pancy ratio R of the system with the rewiring probability 
P = 0.5 is displayed as the system size A'^ is varied. At 
low temperatures (T < 1.5), R first decreases monoton- 
ically with the size A^ and then saturates to the value 
0.5; at high temperatures (T > 2.9), on the contrary, 
R increases monotonically to the saturation value larger 
than 0.5, depending on the temperature. In contrast to 
these monotonic behaviors, at intermediate temperatures 
(1.5 < T < 2.9), the occupancy ratio R behaves non- 
monotonically, exhibiting a maximum at a finite value 
of the system size A^. The height of such a resonance- 
like peak tends to increase as the temperature is raised. 
We have also considered different values of the rewiring 
probability P as well as of the driving frequency fl. It 
is found that as P is increased, the range of the temper- 
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FIG. 1: Behavior of the occupancy ratio R with the system 
size TV at various temperatures T in the absence of long-range 
shortcuts. Lines are merely guides to eyes. 



FIG. 3: The occupancy ratio R vs the system size in the 
presence of the long-range shortcuts (P = 0.5) at various 
temperature T. Lines are merely guides to eyes. 
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FIG. 2: Behavior of the saturation size A'^s with the tem- 
perature T in the absence of shortcuts, exhibiting a linear 
relation between In A^'s and 1/T, with the slope c — 5.6 ± 0.7. 
The dotted line, obtained by the least-square fit, represents 
InNs = 5.6/T + 1.2. 



ature, in which the size-resonance behavior is displayed, 
becomes wider and that the saturation temperature be- 
yond which R shows saturation behavior increases. On 
the other hand, as the driving frequency is increased, the 
position of the resonance peak shifts toward smaller val- 
ues of the system size, thus tending to yield monotonic 
decrease of R. 

Note that this resonance behavior manifests two kinds 
of length scale in the system with long range shortcuts: 
the saturation size Ng and the resonance size at 
which R reaches the maximum. To understand the possi- 
ble relation with the correlation length even in the pres- 
ence of long range shortcuts, we examine the behaviors 
of Ns and Nm, which are displayed in Figs. El 

Figure 0] exhibits the saturation size Ng versus the 
temperature T in the system with P — 0.5. For com- 
parison with the correlation volume described by ~ 
\T — Tc\~'^ 01, the data points are plotted in the log- 



FIG. 4: Behavior of the saturation size Ns with the temper- 
ature T for P — 0.5 at low temperatures (T < Tc), where the 
value Tc — 2.7 obtained from the fitting has been used. The 
least-square fit represented by the dashed line corresponds 
to log TVs = -2.01og(rc - r) -I- 1.6, the slope of which gives 
i? — 2.0. Inset: Behavior of Ns at temperatures higher than 
Tc- The dashed line, obtained from the least-square fit, is 
given by logiVs = -2.01og(r - Tc) + 2.1, where the slope 
again leads to v — 2.0. The error bars have sizes not larger 
than the symbol size. 



arithmic scale, thus fitted to a linear relation between 
IniVg and In jT — Tc\ with the proportionality constant 
(slope) D. From this fitting, we obtain Tc « 2.7, with 
which the slope is estimated to he D = 2.0 ±0.1. Noting 
that in d dimensions, the correlation volume relates with 
the correlation length ^ via £,v ^ S,"^ and the behavior 
^ ~ \T — Tc\^'^ , we thus have D = dv in a d-diniensional 
system. Here it is known that the (effective) dimension 
d oi a. mean-field system should be taken as the upper- 
critical dimension du |14|. leading to v = duV. With 
du — 4: and v = vmf = 1/2 for a mean-field system, 
we conclude that the value v = 2.0 indicates a transition 
of the mean-field nature [3| • We also examine the other 
length scale Nm and show its temperature dependence in 
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FIG. 5: The resonance size Nm together with the saturation 
size Ns vs the temperature T in the presence of shortcuts 
(P = 0.5). The values of Tc, obtained from the fitting, are 
given by 3.1 and 2.7 for Nm and Ns, respectively. The dotted 
line represents the corresponding least-square fit of Nm , which 
is described by logiV„ = -2.01og(rc - T) + 1.2. Thus both 
Nm and Ns result in the same value v — 2.0. 
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FIG. 6: The resonance size Nm vs the temperature T at var- 
ious values of the rewiring probability P. Inset: Behavior of 
Nm with the rewiring probability P at temperature T = 2.2. 



Fig. 1^1 where for convenience iV^, shown already in Fig. 01 
is also plotted. It is shown that both the two length 
scales behave similarly with the temperature, yielding 
essentially the same value of the exponent D — 2.0 ± 0.1 
(see the slopes of the two fitted lines). Accordingly, both 
the two length scales Ns and N„i apparently measure the 
correlation length of the system. On the other hand, the 
fitting parameter Tc for N,n turns out to be 3.1, which 
is somewhat higher than the value 2.7 obtained from fit- 
ting of Ng. At this stage it is difficult to discern unam- 
biguously the difference; more extensive simulations and 
careful analysis are necessary for confirming the origin as 
well as the presence of this discrepancy. 

Finally, we consider systems with different rewiring 
probabilities, and examine how the rewiring probability 
P affects the resonance behavior. FigureElshows the be- 
havior of the resonance size Nm with the temperature at 
various values of P. It is observed that Nm first increases 
slowly with the temperature T then very fast as T ap- 
proaches Tc, which again reminds us of the behavior of 
the correlation length. Behavior of Nm with the rewiring 
probability P at given temperature T — 2.2 is displayed 
in the inset of Fig.|Hl Note the rather fast decrease of Nm 
for small rewiring probabilities (P < 0.5) and the satura- 
tion behavior for large rewiring probabilities {P > 0.5). 
Such saturation behavior has also been reported in the 
synchronization of the system of coupled oscillators on 
a WS network It is also noteworthy that Nm de- 

creases as shortcuts are introduced to the system; the 
shortcuts tend to decrease the optimal system size which 
corresponds to the maximum collective response at given 
temperature and rewiring probability. 

In conclusion, we have investigated the effects of the 
system size on the collective response, measured by 
the occupancy ratio, in the oscillatory field-driven Ising 



model on WS networks. In the purely one-dimensional 
system without long-range interactions, which does not 
undergo a phase transition at any finite temperature, the 
occupancy ratio has been found to display monotonic be- 
havior, not exhibiting a resonance peak. As long-range 
interactions come into the system, on the other hand, 
system-size resonance, characterized by nonmonotonic 
behavior, has been observed to emerge, thus suggest- 
ing a possible relation between the size resonance and a 
finite-temperature phase transition. The resonance size 
at which the occupancy ratio reaches the peak may be re- 
garded as the optimal size of the network, in view of the 
maximum response. It is noteworthy that the optimal 
size as well as the saturation size displays temperature- 
dependent behavior, which is essentially the same as that 
of the correlation length of the system. This suggests the 
interesting possibility of estimating the correlation length 
from the stochastic resonance behavior at various sizes. 
Namely, the size resonance phenomena may be used as 
a tool to measure the correlation length. Note also that 
at given temperature both length scales, the optimal size 
and the saturation size, tend to decrease as the amount 
of long-range interactions is increased. 

As a possible application of the system size resonance, 
we suggest biological systems such as the assembly of 
beta cells which reside in the islets of Langerhans in 
a pancreas |l6|. It is known that the beta cells form 
clusters, each with a finite number of cells, rather than 
gathering together as one unit. Thus speculated is the 
possibility that the function of the beta cells may be op- 
timized via the mechanism of the system-size resonance. 
In addition, the system size resonance behavior may also 
be useful for understanding the formation of the public 
opinion [17i |. 
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